Effects of intermediate bound states in dynamic force spectroscopy 
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We revisit here some aspects of the interpretation of dynamic force spectroscopy experiments. The standard 
theory predicts a typical unbinding force /* linearly proportional to the logarithm of the loading rate r when a 
single energetical banier controls the unbinding process; for a more complex situation of A'^ barriers, it predicts 
at most A'^ linear segments for the /* vs. log(r) curve, each segment characterizing a different barrier. We 
here extend this existing picture using a refined approximation, we provide a more general analytical formula, 
and show that in principle up to N{N + l)/2 segments can show up experimentally. As a consequence the 
interpretation of data can be ambiguous, for the characteristics and even the number of barriers. A further possible 
outcome of a multiple-barrier landscape is a bimodal or multimodal distribution of the unbinding force at a given 
loading rate, a feature recently observed experimentally. 
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I. INTRODUCTION 

The last decades have witnessed a revolution in the methods 
to observe and manipulate single bio-macromolecules or bio- 
molecules complexes. New micromanipulation techniques 
have especially been put forward to probe the folded struc- 
ture of proteins and to quantify the strength of adhesion com- 
plexes (ii 4;i^- An important step in this direction is the 
proposal of the group of Evans to use soft structures to pull 
on adhesion complexes or molecules at various loading rates 
(dynamic force spectroscopy or DFS) (7). Moving the other 
end of the soft structure at constant velocity induces on the 
complex a pulling force that increases linearly in time / ~ rt. 
Measuring the typical unbinding time t* yields an unbinding 
force /* = rt* that depends on the pulling rate r. The typical 
outcome of such experiments is a plot of /* vs. logr com- 
posed by a succession of straight lines with increasing slopes 
(force spectrum). It has been argued and that it is possible 
to deduce the value of some relevant structural parameters of 
the complex by analysing force spectra thanks to an adiabatic 
Kramers picture. This picture consist in considering the un- 
binding process as the thermally activated escape from bound 
states over a succession of barriers along a one-dimensional 
path crossing a mountainous energy landscape ( 7). Within this 
scheme, each straight line of the force spectrum witnesses the 
overcome of an energy barrier and its slope maps the barrier 
to a distance x along the pulling direction. This procedure has 
been shown to yield reasonable values for a few systems, and 
has been conforted by numerical simulations (9). 

Subsequently, theoretical studies have refined the above 
original model, e.g., by inclusion of rebinding events 
study of time dependent loading rates (7), or consideration 
of more complex topographies til3i) and topologies of the en- 
ergetical landscape (Il2i) . 
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In the present paper we explore the potential influence of 
the existence of intermediate bound states on the experimental 
dynamic response of adhesion complexes as probed in DFS. 
To achieve this goal, we first revisit the analysis of the escape 
from a bound state consisting of an arbitrary number of barri- 
ers along a ID path under the application of an external load 
(in line with earlier studies of Strunz et al. (13)), and then dis- 
cuss the implications of this analysis for the interpretation of 
experimental data. In dllH the standard picture is recalled, to- 
gether with its two underlying assumptions. In (IV), we first 
relax the a priori assumption of a deep fundamental bound 
state and provide a general expression that relates the typical 
rupture force to the loading rate (within a single escape rate 
approximation). The practical implications of this new for- 
mula [Eq. ( I15H are discussed, and in particular we comment 
upon intrinsic ambiguities in inferring informnation from a 
[log(r), /*] plot. Then, we show in (V) that in the presence of 
multiple bound states it may be necessary to relax the other as- 
sumption (a single typical rupture force for each loading rate) 
as multimodal rupture force distributions naturally show up, a 
feature recently observed in lipid extraction experiments ( il4l) . 



II. MODEL AND NOTATIONS 

Figure Q] illustrates the energy landscape of a one- 
dimensional escape path with N energy barriers and wells. 
The position and energy of the ith energy well (i = marking 
the fundamental bound state, and 1 < i < — 1 the interme- 
diate ones) are denoted by Xi and Ei, respectively. Similarly, 
the position and energy of the jth energy barrier are denoted 
by Xj and Ej, respectively (where 1 < j < N). For conve- 
nience, without loosing generality, we set xo = and £'o = 
for the fundamental bound state. The unbound "state" is on 
the right hand side of the A^th barrier. If the energy differ- 
ences {Ei ~ Ei) and — Ei) exceed ksT the transition 
rates k~ (and k^) from the ith energy well over the left ith 
barrier (and right (i + l)th barrier, respectively) can be writ- 
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FIG. 1 Sketch of the one-dimensional energy landscape describing 
the unbinding pathway projected along the pulling direction. 

ten according to the Kramers formula 



-(Ei+i-Ei)/kBT 



(1) 
(2) 



A. Single escape rate approximation 

If for any experimentally relevant load / the equilibration 
within the bound states is much faster than the escape to the 
unbound state, the unbinding can be described by a single 
load dependent escape rate fc(/). Following the calculation 
of Evans (T^, if rebinding is negligible (which is the case in 
most experimental situations), the probability P{t) of remain- 
ing in the bound state at time t (the survival probability of the 
bond) then decreases as 



dm 
dt 



-k{rt)P{t). 



(4) 



The solution of this differential equation is P(t) = 
exp[— k{rt') dt']. The probability density for unbinding 
between times t and t + At is Pt{t) = —dP{t)/dt = 
k{rt)P{t), from which, after changing the variable from t to 
/, one gets the probability density for the distribution of the 
unbinding force: — {l/r)k{f)P{f /r). The typical un- 

binding force /* is defined as the peak of this probability den- 
sity: dpf{f)/df\f=f* = 0, which yields the simple formula 



where loq is a typical attempt frequency, a,; and aj 
metric factors characterizing the shape of the ith energy well 
and jth energy barrier, respectively. Note that there is no tran- 
sition from the fundamental bound state to the left, therefore, 
k, ^ 0. 

We assume throughout the paper that the energy wells and 
barriers are sharp, so that for any loading force / their loca- 
tions remain constant, and their energies change as Ei{f) = 
Ei{Q) - fx, and Ej{f) = Ej{Q) ~ /%. To simpHfy the no- 
tations, wherever the argument of the energies and transition 
rates is omitted, a loading force / is implicitly assumed. 

Finally it will prove convenient to introduce a few com- 
pact notations. For any < i < j < we denote the 
distance between the ith well and the jth barrier (on the 
right) by Axij — xj — Xi, and their energy difference by 
AEij = Ej — Ei. We also define a formal (effective) rate 
constant from the ith well over the jth barrier on the right as 



(3) 



Obviously AEij and kij implicitly depend on /, whereas 
Axij are constants given the assumption of the previous para- 
graph. 



dr(/) 



d/ 



/=/* 



(5) 



where r(/) = l/k{f) denotes the load dependent mean es- 
cape time. This formula gives the loading rate r at which the 
typical unbinding force is /*. For practical purposes it is of- 
ten necessary to invert this relation to, e.g., predict the typical 
unbinding force for an experimentally imposed loading rate. 

To set a reference for further comparison, we explicitely 
invert the above relation in case of a single barrier, i.e. when 
t(/) — t(0) exp(— /xi), and obtain 



Xl 



rT(0)xi 



kBT 



(6) 



As mentioned in the introduction, the escape over a single bar- 
rier results in a single straight line in the force spectrum. The 
experimental observation of a linear segment consequently 
gives hints as to the structure of the energy landscape, in par- 
ticular the slope of the segment permits to deduce a distance 
Xl between the energy well and the barrier. 



B. Deeply bound fundamental state approximation 



III. STOCHASTIC KINETICS OF UNBINDING UNDER 
EXTERNAL FORCES: STANDARD DESCRIPTION AND 
CORRESPONDING APPROXIMATIONS 

We first recall the standard description of the "force spec- 
trum", which relies on two major assumptions, namely the 
single escape rate (SER) and the deeply bound fundamental 
state (DBFS) approximations. 



Assuming further that the fundamental bound state is much 
deeper than the intermediate ones: Ei{f) — Eo{f) ^ fcerfor 
any experimentally relevant load / (i.e. before unbinding has 
statistically almost certainly occurred, see Fig.|2a.), Evans has 
shown that the mean escape time from the fundamental bound 
state to the unbound state is well approximated by 



N 



if) 



-fAxoj/ksT 

E^^-^^ (7) 
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FIG. 2 Sketch of two energy landscapes with one intermediate well. 
Dotted drawings: no external force. When a constant force is applied 
energies are lowered by fx (dashed lines), the resulting landscapes 
appear in solid lines. The dotted arrows indicate which pair of wells 
and barriers control the kinetics at zero load. The solid arrows indi- 
cate the new limiting effective escape process at higher forces, (a) 
The escape from the fundamental bound state remains the limiting 
process whatever the pulling force, (b) The escape form the interme- 
diate bound state energy becomes the limiting process at high forces. 

This allows one to obtain an explicit r vs. /* relationship by 
plugging Eq. Q into Eq. (|5}, which yields the compact for- 
mula 



N 

E 



-f'Axo,,/ksT 



(8) 



This equation predicts a spectrum /* vs. log(r) consisting of 
a succession of at most N segments of increasing slopes, each 
of which yielding an information Axqj about an intermediate 
barrier. 



IV. BEYOND THE DEEPLY BOUND FUNDAMENTAL STATE 
APPROXIMATION 

In this section we relax the approximation made in the last 
subsection (IIIB), generalize accordingly equations (7) and 
(9), and discuss the experimental implications of this gener- 
alization. 



A. Refined theory 

In general, it is possible that for large enough forces one 
or more of the intermediate bound states become deeper than 
the fundamental bound state before unbinding has occured 
(see Fig. 12 b.). In such cases the above DBFS approxima- 
tion breaks down. However, we show below that it is still 
possible to compute rather simply the escape time r from the 
"bound state" to the "unbound state", provided we maintain 
the assumption of a single escape rate l/r(/). 

Let us put the system into its fundamental bound state, and 
let it evolve according to the transition rates given in Eqs. Q 



and (|2jl. Whenever the system gets into the unbound state (by 
making a transition over the outermost barrier) let us place it 
back into the fundamental bound state. The stationary state 
of an ensemble of such systems is characterized by a proba- 
bility current, which is constant everywhere and equal to 1/r 
by definition. To calculate t we have to solve the following 
system of equations: 



P,k+ ~ P,+iK^, = l/r 



N-l 



0<i<N-2, (9) 
(10) 

(11) 



where Pi denotes the probability of being in the ith bound 
state (0 < i < — 1). The first N equations describe the 
probability current over each of the N barriers, and the last 
equation is just the normalization condition. These + 1 
linear equations uniquely determine the + 1 variables (Pi 
and r), and can be solved easily in a recursive way. First, 
Pn-it can be expressed from Eq. ( II 0> . and then Pjv_2T, 
Pqt recursively from Eq. (|9} yielding 



1 

N 



k 



■"-Af-l 



F — 



(12) 



j=i+i 



where Eqs. Q, (|2j, and the definition Q have been used. 
Note that because the ki^j are only formal definitions, con- 
structed as products and ratios of the single-barrier rates ([TJ 
and (|2ji, they are meaningful even if AEij < 0. From the 
normaUzation (^J one can easily express t as 



N-1 N 



E^^-E E 



1=0 



1=0 



ki 



(13) 



The sum is dominated by the smallest effective rates, which 
are the bottlenecks of the unbinding process. Consequently, 
this formula remains a good approximation for r even if some 
of the barriers disappear at big loads, because the correspond- 
ing formal transition rates make negligible contributions. By 
indicating the load force / explicitly, we arrive at 



rif) 



N 

E 



1 



=0 j"=i+l *'-'^*' ' 



N~l N 

E E ■ 

i=0 j=i+l 



(14) 

which generalizes (7). An analytic formula can be given for 
the /* vs. r relationship by plugging Eq. MA\ into Eq. (|5} : 



N 



E 



Ax 



2J 



^ ^ k^T 

i=0 j=i+l 



-f-Axij/ksT 

~k~{0) 



(15) 



This generalization of equation (8) is one of the main results 
of this paper. Let us briefly comment on immediate features 
of this new formula. 
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First, Eq. (|8} is easily recovered from (15) assuming a 
DBFS. Indeed, the assumption Ei{f) 3> EQ{f) implies 
ko.j 3> ki,j, if i > [see Eq. (|3}] and therefore, the relation 
Pi / Po <C 1, if i > is deduced from Eq. (I12> . The probabil- 
ity to find the system in the fundamental bound state is close 
to 1. So, the sum over i in Eq. MA\ is dominated by the con- 
tributions of the effective escape rates from the 0th well only. 
Finally, the sum over i (labeling the intermediate states) is re- 
duced to its sole first term too, and Eq. M5\ becomes identical 
to Eq (|H). 

Second, each of the N{N + l)/2 terms of Eq. M5\ alone 
would yield a straight line in the /* vs. log(r') plot. However, 
at any loading rate the highest force value (the uppermost seg- 
ment, corresponding to the most difficult transition) limits the 
unbinding process, therefore, the /*(?') curve is expected to 
closely follow the upper envelope of these segments [see Fig. 
13 a)]. Depending on the position of the lines, this upper enve- 
lope can consist of up to N{N + 1) /2 linear segments. 

Third, this last point is clearly at odds with the prediction 
within the DBFS approximation. Indeed assuming a DBFS 
corresponds to forbidding the display in the force spectrum of 
the N{N — l)/2 segments corresponding to the probing of the 
escape from an intermediate bound state (see Fig.|2] b.). 



B. Practical implications: Ambiguity In the determination 
of "structural" parameters 

We now insist on some practical implications of the above 
general description. We do not attempt a full inspection of 
all the possible dynamic responses of arbitrarily complex sys- 
tems, but rather focus on two simple examples in order to 
stress that the main features of the energy landscape can in 
general nofbe unambiguously inferred from [log(r), /*] plots. 
To emphasize the experimental relevance of this discussion, 
we use for the parameters values comparable to those ob- 
served in experimental systems. Specifically, we take the ge- 
ometric factors Q!iS and SjS to be all equal to 1, cjo = 10* s^^ 
and fceT = 4 x IQ-^ij. 



1 . Ambiguity in determining tlie barriers positions 

Fig|3](a) and Fig|3](b) display two force spectra as obtained 
from Eq. (I15> . Both correspond to energy landscapes with 
two barriers. Though the two [log(r),/*] plots are almost 
identical they are related to very different sets of values for 
the energy levels and positions (along the pulling direction) of 
the wells and the barriers. 

Fig|3](a) corresponds to the situation where the standard 
picture to account for the two segments is well suited (7). 
At low force, the escape from the fundamental 0th state over 
the outermost barrier is the limiting process. The slope of 
the first segment is proportional to ksT /x2- For the highest 
forces (above ~ 30pN) the energy of the external barrier is re- 
duced below El and the deepest bound state remains located 
at xq = 0. The process that mostly impedes the unbinding is 
the overcome of the innermost barrier i?i with arate fco,i. The 
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FIG. 3 Two very similar force spectra corresponding to different en- 
ergy landscapes with one intermediate well. 

Curves in solid lines: Force spectra plotted using Eq. <15> . It closely 
follows the upper envelope of the straight lines corresponding to the 
transitions: fco,2 (dotted line), fco.i (dashed line) and k\^2 (dash- 
dotted line). 

Insets: shape of the energy landscapes at low and high forces, each 
drawn with the same line style as the straight line associated with the 
limiting transition. 

Parameter values: (a) (3;i,i5i) = (Inm, llfcBr), (x\,E\) = 
(1.5nm,8fcBr), and (s2,i?2) = (2 nm, 20 ^7), (b) Si) = 
(0.5nm, 12fcBr), {xi.Ei) = (lnm,9fcBr), and {x2,E2) = 
(2nm,20 /cbT). 

slope of the curve is now larger and proportional to fesT/xi. 

Fig. |3](b) corresponds to an energy landscape for which the 
above explanation is inappropriate. At low force the unbind- 
ing kinetic is controlled by the escape from the fundamental 
state over the outermost barrier again. But, for pulling forces 
larger than ^ 30pN this outer barrier remains the highest [see 
inset in Fig|3](b)]. However the slope of the spectrum in- 
creases as in the (a) case. The reason is that the deepest (and 
most occupied) bound state is now located ?A x — xi and 
the presence of the second segment actually witnesses the es- 
cape from this intermediate state to the unbound state with a 
rate fci_2. The value of the second slope scales therefore with 
ksT/ {x2 — xi). Since the escape rate ki^2 in the (b) case is 
equal to the escape rate fco.i in the (a) case, the two spectra in 
Fig-EJturn out to be indistinguishable and cannot be a priori 
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FIG. 4 Two very similar force spectra corresponding to energy land- 
scapes with different numbers of intermediate wells. 
The rule of the line styles is the same as in Fig.|3| (a) dotted lines: 
fco,3, dash-dotted lines: fco,2, dashed lines: fco,i- (b) dotted lines: 
fco,2, dash-dotted lines: k\.2, dashed lines: fco,i. 
Parameter values: (a) {x\^ E{) = (0.6 nm, 14 fcer), (a;i,i5i) — 
(0.7nm, 12fcBr), (s2,i?2) = (1.1 nm, 19 /cbT), (a;2,-E;2) = 
(2nm, lefcaT) and (£3,^3) = (2.5 nm, 24 fcBT). (b) (Si,^i) = 
(l.lnm, igfeeT), (si.^i) = (1.9 nm, 10 fcaT), and (^2, S2) = 
(2.5 nm, 24 fceT). 

associated with one of the two possible landscapes. 



2. Ambiguity in determining tlie number of barriers 

After having shown with the simple example above that am- 
biguity can exist in determining distances from dynamic force 
spectra, we show here that even more strikingly it is impos- 
sible in general to assess the number of wells and barriers. 
Again we use a simple example to do so. 

Fig. |4]displays two force spectra obtained using Eq. ( I15t . 
They are both well approximated by a succession of three 
segments with increasing slopes. Again, the two [log(r), /*] 
curves are very similar although they are constructed from 
landscapes that do not even comprise the same number of 
peaks and wells. 

In Fig. 13(a) the three segments describe the escape from the 
same fundamental state over the three distinct energy barriers. 



The larger the pulling force the closer the limiting barrier to 
the fundamental state [see inset in Fig.0](a)]. 

In Fig. 13(b), the landscape consists of only two barriers. 
However, the force spectrum reveals that three different es- 
cape processes can limit the unbinding kinetic. At low forces 
(/ ^ 50pN) the two observed Unear segments results from 
the escape form the fundamental state over the two peaks at 
x\ and X2 respectively. Conversely, at high forces it is the es- 
cape from the deeply lowered intermediate state over the outer 
baiTier that determines the escape rate (see inset in Fig. |3(b), 
drawing with dash-dot line). With the chosen parameters the 
effective rate fco,3^ ^0,2 and fco,i in Fig- 4(a) case coiTespond 
respectively to fco,2, ^0,1 and fci.2 in Fig. 4(b) case. Thus the 
two plots are indistinguishable and cannot be used to predict 
the number of barriers along the ID escape path. 

In conclusion of this subsection, we suggest great care in 
infeiTing features of the underlying energy landscape from dy- 
namic force spectroscopy experiment. Our generalized equa- 
tion may be helpful in dealing with the coiTesponding ambi- 
guity as it allows (with some work) to generate various land- 
scapes that can account for the observed data, whereas Eq. (|8jl 
can only yield a single set of parameters ( e.g. those used for 
the plots in Fig.|3l(a) and Fig.|3(a)). 



V. BEYOND THE SINGLE ESCAPE RATE 
APPROXIMATION: MULTIMODAL UNBINDING FORCE 
DISTRIBUTIONS 

Up to this point we have been considering a generalized 
theory in which the deeply bound fundamental state (DBFS) 
approximation is dropped, but the unbinding is still approxi- 
mated as a simple first-order escape process. Indeed, the va- 
lidity of Eq. dlSI l relies on the assumption that at any moment 
the distribution of the populations of the bound states can be 
well approximated by the distribution corresponding to a ho- 
mogeneous stationary cuiTent. 

This is, however, not always the case. As we stated earlier, 
the sum of the l/kij{f) terms in Eq. il4\ is dominated by 
the smallest effective rate constant ki'j'{f) coiTesponding to 
the slowest effective transition. A consequence of this is that 
all the bound states located to the left of barrier j' are close 
to equilibrium (because of the slow outflow over barrier j'), 
and the population of any state located to the right is negli- 
gible (because they practically belong to the unbound state). 
Now, If the slowest transition rate changes from to 
if) as the loading force / is increased, and if j" < j', 
then a considerable population might remain in the intermedi- 
ate bound states between the new and the old limiting barriers, 
j" and j' respectively. This residual population is incompat- 
ible with the new stationary current dominated by kin ^jii{f), 
and must escape in a different way, yielding a secondary max- 
imum of the unbinding force distribution (see Fig.lSji). 

The escape of the majority of the population (located to the 
left of the new limiting baiTier j") can still be characterized 
by Eq. ( I14> of our generalized theory. On the other hand, we 
have to slightly modify this formula to describe the escape of 
the residual population (trapped between the new and old lim- 
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FIG. 5 A scenario that yields a multimodal unbinding force distri- 
bution.(a) Three snapshots of the energy landscape for the pulling 
forces: / = OpN (dotted line), / = 18 pN (dahs-dotted line) 
and / = 60 pN (dashed line). The three arrows indicate the cor- 
responding most difficult transitions. Parameters values: (xi, Ei) — 

(0.5nm, lefeBT), ixi,Ei) = (1.1 nm, 4 fcBT), and (^2, i?2) = 
(2.3 nm, 24 ksT). (b) Force spectrum (full line) associated with the 
landscape described in (a) plotted using Eq. <15> . The full probabil- 
ity density P/(/) for unbinding at force / is plotted in gray scale in 
the inset. p/(/) has been obtained following the procedure described 
in IT^ . The circles in the main plot represent the local maxima of 
the unbinding force distribution Pf{f) 

iting barriers). Since j" is the limiting barrier now, almost the 
entire residual population can escape without ever jumping 
backward over barrier j". Therefore, for the residual popula- 
tion we can consider barrier j" as a reflecting boundary, and 
describe the escape by our general theory in this modified po- 
tential. Eq. J14l i. e.g., changes accordingly: 



N-l N 



-if) = E E 



1 



N-l N 

E E 

i.=j" j=«+i 



(16) 



Consequently, the absolute maximum of the unbinding 
force distribution always follows the upper envelope of the 
N{N + l)/2 lines, however, some secondary maxima might 
also appear at lower forces, which follow the upper envelope 
of only a subset of the lines [comprising {N — j"){N — j" + 
l)/2 elements]. Such secondary maxima of a multimodal 
force distribution give important information on the internal 
structure of the energy landscape of the unbinding path, and 
makes the determination of the number and positions of the 
energy wells and barriers less ambiguous. It is actually a 
nice achievement of our generalized theory to be able to make 
sense of the segments of secondary maxima in a unique frame 
for fitting parameters (see e.g. Fig. 5, where the segment 
corresponding to the secondary maximum corresponds to the 
transition from the well 1 over the barrier 2, a step neglected 
in the DBFS approximation). The possibility of a bimodal dis- 
tribution for the case of a two-state system has already been 
reported by Strunz et al. and our description system- 

atizes and generalizes their findings. 

To provide a simple illustration for the somewhat formal 
discussion above, we also focus on a system consisting of two 
bound states, as depicted in Fig.|5](a). Increasing the force, 
the limiting transition rate changes from fco.2 to fci 2 and then 
to fco,i- In the range of the loading rate r between about 10^ 
and 10^ pN/s the intermediate bound state (1) has enough time 
to accumulate a large population, which is then flushed by the 
/ci_2 transition before the transition fco.i flushes the rest from 
the fundamental bound state (0). In the range above 10^ pN/s 
the intermediate bound state (1) cannot accumulate much of 
the population, but it still possesses a small fraction of the 
initial equilibrium distribution, which is again flushed by the 
/ci 2 transition first. 

Very recently the group of Evans actually reported the ex- 
perimental occurrence of a bimodal force distribution 
The corresponding experiment consisted in pulling on "diC14 
PE" lipids from a bilayer made of "C18:0/l PC" lipids. With 
the help of our generalized theory, Evans and Williams were 
able to fit their data and interpret the results in terms of an 
energy landscape with two barriers [personal communication, 
seeRef. [3] in Ref. Q]. 



VI. CONCLUSION 

In this paper, we have revisited the standard theory used 
to account for the dynamic response of molecular stickers. 
Our refined description, valid for an arbitrarily complex one- 
dimensional energy landscape, has allowed us to highlight 
several practical consequences of the diversity of the possible 
unbinding scenarios. For example several markedly different 
energy landscapes can yield the same rupture force distribu- 
tion. To resolve this ambiguity other experimental technics, 
e.g. flow chamber experiments f\3), are then required. We 
have also identified the physical origin of multimodal unbind- 
ing force distributions and shown how their analysis provides 
informations on the unbinding pathways. 
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